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Abstract The evolution of ozone is examined in the latest version of the Goddard Earth Observing System 
Chemistry-Climate Model (GEOSCCM) using old and new ozone-depleting substances (ODS) scenarios. This 
version of GEOSCCM includes a representation of the quasi-biennial oscillation, a more realistic implementation 
of ozone chemistry at high solar zenith angles, an improved air/sea roughness parameterization, and an extra 
5 parts per trillion of CH 3 Br to account for brominated very short-lived substances. Together these additions 
improve the representation of ozone compared to observations. This improved version of GEOSCCM was used 
to simulate the ozone evolution for the A1 201 0 and the new Stratosphere-troposphere Processes and their Role 
in Climate (SPARC) 2013 ODS scenario derived using the SPARC Lifetimes Report 2013. This new ODS scenario 
results in a maximum Cl tot increase of 65 parts per trillion by volume (pptv), decreasing slightly to 60 pptv by 
21 00. Approximately 72% of the increase is due to the longer lifetime of CFC-1 1 . The quasi-global (60°S-60°N) 
total column ozone difference is relatively small and less than 1 Dobson unit on average and consistent with the 
3-4% larger 2050-2080 average Cl y in the new SPARC 201 3 scenario. Over high latitudes, this small change in Cl y 
compared to the relatively large natural variability makes it not possible to discern a significant impact on ozone 
in the second half of the 21st century in a single set of simulations. 


1. Introduction 


Ozone is a critically important trace gas in the Earth climate system. It is closely connected to the dynamics of 
the atmosphere as well as the climate [Gillett and Thompson, 2003; McFarlane, 2008]. Since 1994 the World 
Meteorological Organization (WMO) has conducted a careful review of observations of ozone and related 
constituents and state-of-the-art simulations every 4 years, most recently producing Scientific Assessment of 
Ozone Depletion : 2010 [World Meteorological Organization (WMO), 201 1]. While many factors influence the 
evolution of ozone, changes in the levels of ozone-depleting substances (ODS) and the Brewer-Dobson 
circulation are important drivers of ozone change [Oman et al., 2010]. 

Recently, a new assessment of the lifetimes of many ODS was undertaken by the World Climate Research 
Program's project Stratosphere-troposphere Processes and their Role in Climate (SPARC) [SPARC, 2013], 
leading to a new ODS scenario based on the latest scientific understanding. One of the larger changes was to 
the lifetime of CFCI 3 (CFC-1 1, F1 1 hereafter), which was reassessed to be 52 years from 45 years based on Kaye 
et al. [1994]. This longer lifetime results in higher atmospheric loading of F1 1 in the 21st century. 1 /elders and 
Daniel [2014] produced a new mixing ratio-based boundary condition ODS scenario based on the reassessed 
F1 1 lifetime and a number of other species. They found a higher level of equivalent effective stratospheric 
chlorine (EESC), which is the total effective amount of chlorine and bromine in the stratosphere [Daniel et al., 
1995; Newman etal., 2007], mainly during the second half of the 21st century as compared to the previous Al 
2010 scenario [WMO, 2011], 

In addition to halogen loading, possible Brewer-Dobson circulation changes in response to climate change 
[Butchart et al., 201 1] could have important impacts on ozone and other atmospheric constituents [Oman et al., 
2010]. The evidence for changes in the Brewer-Dobson circulation has so far been mixed, with some studies 
showing little if any change [Engel et al., 2009; Ray et al., 201 0] to others suggesting significant tropical lower 
stratospheric changes [Randel and Thompson, 2011; Kawatani and Hamilton, 201 3] over the most recent decades. 

The quasi-biennial oscillation (QBO) is a leading mode of variability in the tropical stratosphere [Baldwin etal., 
2001], The QBO modulates deep tropical upwelling and temperature, which strongly impacts the variability 
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Table 1 . Simulations Used in This Study 


Name 

Years 

SST/SIC 

GHG 

ODS 

+5 Br y 

QBO 

Figures 

GEOS-4_o 

1 960-2004 

Obs. a 

Obs. a 

Obs a 

No 

No 

1 and 2 

GEOS-5 old o 

1 960-2005 

Obs. a 

Obs. a 

Obs a 

No 

No 

1 and 2 

GEOS-5 new o 

1960-2012 

Obs. a 

Obs. a 

Obs a 

Yes 

Yes 

1 and 2 

GEOS-4_m 

1 972-2098 

CCSM3 b 

A1B d 

Al 2002 f 

No 

No 

3 and 4 

GEOS-5 old m 

1 960-2099 

CCSM3 b 

A1B d 

Al 2006 9 

No 

No 

3 and 4 

GEOS-5 new m 

1 960-2099 

CESM1 c 

RCP6.0 e 

SPARC 201 3 h 

Yes 

Yes 

3 and 4 

P Al 2010 

1960-2012 

Obs. a 

Obs a 

Obs a 

Yes 

Yes 

6, 7, and 8 

PF Al 2010 

1 960-2099 

CESM1 c 

RCP6.0 e 

Al 2010' 

Yes 

Yes 

5b, 6, 7, 8, 9, and 10 

PF_SPARC_201 3 

1 960-2099 

CESM1 c 

RCP6.0 e 

SPARC 201 3 h 

Yes 

Yes 

5b, 6, 7, 8, 9, and 10 


a Obs. refers to Observations. 

“Community Climate System Model version 3 (CCSM3). 
c Community Earth System Model version 1 (CESM1). 

d A1 B refers to the Intergovernmental Panel on Climate Change (IPCC)/Special Report Emission Scenario A1 B. 
e RCP6.0 refers to the Intergovernmental Panel on Climate Change/Representative Concentration Pathways Scenario 
6.0 (IPCC/RCP6.0). 

f A1 2002 refers to the WMO 2002 Ozone Assessment ODS scenario A1. 

9 A1 2006 refers to the WMO 2006 Ozone Assessment ODS scenario A1. 

h SPARC 2013 refers to the new ODS scenario derived using the SPARC lifetimes Report 2013. 

1 A1 2010 refers to the WMO 2010 Ozone Assessment ODS scenario A1. 


of ozone and many other constituents and influences the mean state of ozone as well as other constituents 
[Hurwitz et al., 2013]. Simulation of a realistic QBO is a new feature for the Goddard Earth Observing System 
version 5 (GEOS-5) general circulation model and when coupled with a chemistry module allows the 
examination of the QBO impact on trace gas constituents [Hurwitz et al., 2013]. This and other improved 
capabilities in the latest version of GEOS Chemistry-Climate Model (GEOSCCM) are examined here, focusing 
mainly on total column ozone. 

The main goals of this study are (1) to document the improvements in ozone in GEOSCCM, (2) to demonstrate 
the changes due to improvements in GEOSCCM, and (3) to compare the evolution of ozone in the new ODS 
scenario to that obtained with the previous ODS scenario. Section 2 describes the improved version of 
GEOSCCM, the boundary condition scenario, and observations of ozone that are compared with simulations. 
Section 3 discusses the improvements to ozone in GEOSCCM and compares the evolution of ozone in the old 
and new ODS scenarios. Section 4 summarizes the main results and gives concluding remarks. 

2. Model Description, Scenarios, and Measurements 

We discuss improvements to the simulation of ozone and examine the evolution of ozone to new and old 
ODS scenarios in the GEOSCCM. This latest version of GEOSCCM couples the GEOS-5 [Rienecker et al., 2008; 
Molod et al., 201 2] general circulation model to a comprehensive stratospheric chemistry module. Version 4 
of GEOS (GEOS-4) with stratospheric chemistry was discussed in detail in Pawson et al. [2008], and any 
changes that significantly impact ozone and related species since that paper are detailed below. The GEOS-5 
simulations described in this study all use a horizontal resolution of 2° latitude and 2.5° longitude with 72 
vertical layers up to 0.01 hPa (80 km); the GEOS-4 simulation has the same horizontal resolution and vertical 
domain but only 55 layers. Table 1 gives an overview of the simulations considered in this study, including 
simulation length, forcings applied, whether or not there was a QBO, and the figures associated with 
each simulation. 

GEOSCCM has been extensively evaluated using process-oriented diagnostics in both Chemistry-Climate 
Model Validation 1 (CCMVal-1) [ Eyring et al., 2006] and CCMVal-2 [SPARC Chemistry-Climate Model Validation 
(i CCMVal ), 2010]. GEOSCCM has performed well in both chemical- and transport-related processes [SPARC 
CCMVal, 2010; Strahan et al., 201 1] Douglass et al., 2012], but noted deficiencies include the following: the lack 
of a QBO, late breakup of the Antarctic vortex, and high bias in high-latitude ozone. The latest version of 
GEOSCCM includes significant improvements to most of these deficiencies. 

GEOSCCM now internally generates the QBO [Hurwitz et al., 2013] with realistic periodicity and magnitude as 
a result of improvements to the nonorographic gravity wave drag scheme. Revisions to the air/sea roughness 
parameterization [Garfinkel et al., 201 1, 2013] improved the surface wind speeds, especially over the Southern 
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Month 



c) Seasonal Average (50°N-55°N) 1995-2004 



Month 

Figure 1 . The annual cycle of total column ozone averaged between 1995 and 2004 over (a) 75-80°S, (b) 20°S-20°N, and (c) 
50-55°N for GEOS-4_o (green curve), GEOS-5_old_o (red curve), GEOS-5_new_o (light blue curve), and both ground-based 
(solid black curve) and satellite-derived (dashed black curves) measurements. The gray-shaded areas around the ground- 
based data and color whiskers on the simulations are ±2 standard deviations. 
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Ocean, and reduced the late breakup of 
the Antarctic polar vortex. Significant 
improvements in the simulation of ozone 
have also resulted from allowing the 
transition from day/night chemistry to 
occur at a maximum solar zenith angle of 
94° rather than the 90° used previously. An 
extra 5 parts per trillion (ppt) of CH 3 Br is 
now included to represent very short-lived 
brominated substances and is in better 
agreement with observations [Liang et a!., 
2010]. Observed/modeled stratospheric 
sulfate surface area densities [Eyring et al., 
201 3] are used in the new simulations from 
1960 to 201 1 with 201 1 values used for 
years after 2011. Previous simulations 
assumed an unperturbed stratospheric 
aerosol background state from 1979. 

The old (Al 2010) and new (SPARC_2013) 
ODS scenarios are described in Velders and 
Daniel [2014] and are both constrained by 
observations in the past (up to -2008). All 
of the new GEOSCCM simulations include 
greenhouse gas (GHG) concentrations 
from the Representative Concentration 
Pathway (RCP) 6.0 indicating 6.0 W/m 2 
radiative forcing by 2100 [ Meinshausen 
et al., 2011; Moss et al., 2010]. While the 
RCPs also provide scenarios of ODS, we are 
not using them in this study. Sea surface 
temperature (SST) and sea ice 
concentrations (SIC) were prescribed from 
an AR5 simulation using the Community 
Earth System Model version 1 (CESM1) 
from 1 960 to 2099 [Cent et al., 2011] forced 
with the same RCP6.0 GHG scenario. 

Ozone observations used in this study come from satellites and ground-based observing networks. The 
satellite observations are from the Total and Profile Merged Ozone Data Set version 8.6 [Bhartia et al., 2013; 
Stolarski and Frith, 2006; Ziemke et al., 2005] (http://acbd-ext.gsfc.nasa.gov/Data_services/merged/index. 
html), and the ground-based total column ozone (TCO) measurements are updated from Fioletov et al. 
[2002, 2008], 

3. Results and Discussion 

3.1. Ozone and Tropical Upwelling in GEOSCCM 

We first focus on how the ozone and tropical upwelling have evolved between previous versions of 
GEOSCCM and the current version. For Figures 1-4, the green curve (labeled GEOS-4_o or GEOS-4_m) is from 
the version used in CCMVal-1, red curve (GEOS-5_old_o or GEOS-5_old_m) is from CCMVal-2, and the light 
blue curve (GEOS-5_new_o or GEOS-5_new_m) is the latest version that includes the QBO and additional 
changes discussed in section 2. GEOS-5_new is the same version that is used to produce the two ODS 
scenarios simulations discussed later. 

First, we examine the seasonal cycle in the Northern Hemisphere (NH) midlatitudes, tropics, and Southern 
Hemisphere (SH) high latitudes compared to both ground- and satellite-based estimates of TCO. We use 


a)63°s-90°s October 




Figure 2. Total column ozone (DU) over area-weighted polar cap 
averages for (a) 63-90°S during October and (b) 63-90°N during 
March for GEOS-4_o (green curve), GEOS-5_old_o (red curve), GEOS- 
5_new_o (light blue curve), and both ground-based (solid black curve) 
and satellite-derived (dashed black curves) measurements. 
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Annuol Avg. 20°S-20°N w* (mm/s), 70 hPa 



Time (years) 

Figure 3. Annual average tropical (20°S-20°N) residual vertical velocity 
(mm/s) at 70 hPa for GEOS-4_m (green curve), GEOS-5_old_m (red 
curve), and GEOS-5_new_m (light blue curve). GEOS-4 simulation is 
from 1 972 to 2098, and the GEOS-5 simulations are from 1 960 to 2099. 
The thick lines are output smoothed to remove interannual variability. 


thel 995-2004 time period for comparison 
since it is the common period among all 
simulations and includes a well-developed 
ozone hole. Figure la shows that the TCO 
over 75-80°S was very similar for GEOS-4_o 
and GEOS-5_old_o with both biased high 
compared to observations for all months 
other than October-December. The bias 
was particularly large during the ozone 
hole onset in August and September. A 
significant portion of this high bias has 
been reduced in the GEOS-5_new_o 
version from two particular changes, 
allowing the ozone chemistry to continue 
at high solar zenith angles and the 
inclusion of an extra 5 ppt of Br y . The extra 
Br y had its largest impact during August 
and September, whereas the high solar 
zenith angle change reduced ozone year 

round. Lamago et al. [2003] noted a similar reduction in ozone when the high solar zenith angles were allowed 
to impact chemistry in the European Centre/Flamburg 4 model. Figure 1b shows that the tropical (20°S-20°N) 
ozone is well represented in all versions of GEOSCCM with a small decrease in ozone between GEOS-5_old_o 
and GEOS-5_new_o. The NFi midlatitudes (50-55°N) are shown in Figure 1c with a general high bias outside of 
boreal summer. This high bias is partially reduced in GEOS-5_new_o compared to GEOS-5_old_o but similar to 
the GEOS-4_o version. 

Next, we compare the past changes in polar total column ozone for October and March for the SH and NFI, 
respectively. Figure 2a shows the area-weighted polar cap (63-90°S) average TCO for October in observations 
(black solid and dashed curves) and the three versions of GEOSCCM. The main difference in simulations can 
be seen before 1980 with some reduction in the high bias in the GEOS-5_new_o run. A similar but smaller 
reduction in high bias is also apparent in the NFI polar region (63-90°N) in March; however, some biases still 
remain. In the SFI polar region the root-mean-square (RMS) difference between 1964 and 2004 is 67 Dobson 
units (DU) for GEOS-4_o and GEOS-5_old_o and is reduced to 39 DU for GEOS-5_new_o compared to ground- 

based TCO measurements. In the NFI polar 
region the RMS difference is 53 DU for 
GEOS-4_o and 62 DU for GEOS-5_old_o 
and is reduced to 44 DU for GEOS- 
5_new_o. None of these simulations 
reproduces the string of low-ozone years 
seen in observations in the mid-1990s. 

The residual vertical velocity (w*) from 20°S 
to 20°N at 70 hPa can be used as a proxy 
for Brewer-Dobson circulation changes. 
Tropical upwelling changes affect the 
distributions of ozone [Li et al., 2009; Oman 
et al., 201 0] and other trace gases so that 
any change in the predicted rate of 
tropical upwelling increase can impact 
their distributions. This next set of 
simulations is with the same three versions 
of GEOSCCM considered in the previous 
figures but with modeled SST and SIC. In this 
study GEOSCCM is not coupled to an 
interactive ocean, so SST and SIC are 



Figure 4. Annual average quasi-global (60°S-60°N) total column ozone 
(DU) for GEOS-4_m (green curve), GEOS-5_old_m (red curve), GEOS- 
5„new_m (light blue curve), ground-based measurements (black solid 
curve), and satellite-based measurements (black dashed curve). GEOS- 
4_m simulation is from 1972 to 2098, the GEOS-5 simulations are from 
1 960 to 2099, and the ground-based observations are from 1 964 to 201 1 . 
The thick lines are output smoothed to remove interannual variability. 
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prescribed from Community Climate 
System Model version 3 (CCSM3) and 
CESM1 coupled atmosphere ocean 
simulations (see Table 1). The GEOS-5 
simulations span 1960-2099, while the 
GEOS-4_m simulation was integrated from 
1 972 to 2098. While all simulations show an 
increase in tropical upwelling over time 
(Figure 3), the GEOS-4_m simulation shows 
the largest increase. This increase is 
reduced in subsequent versions of the 
GEOSCCM. It is important to note that while 
the GEOS-4_m and GEOS-5_old_m 
simulations use the same SST and A1 B GFIG 
scenario, GEOS-5_new_m uses SSTs and 
GFIGs consistent with the RCP6.0 scenario. 
The differences in tropical SSTs between 
GEOS-5_old_m and GEOS-5_new_m (not 
shown) are consistent with the differences 
in tropical upwelling in these two 
simulations. There is also little apparent 
difference in the upwelling interannual 
variability at 70 hPa between GEOS- 
5_old_m (no QBO) and GEOS-5_new_m 
(QBO) over to 20°S-20°N region. This point 
will be examined in more detail later 
comparing ozone changes in the deep 
tropics to the broader tropics. 

The annual average quasi-global (60°S- 
60°N) TCO for the same simulations is 
shown in Figure 4. GEOS-4_m and GEOS- 
5_old_m produce values that are high 
biased compared to ground- and satellite-based observation of the recent past. The high bias is greatly 
reduced in GEOS-5_new_m. GEOS-5_new_m includes observed/modeled sulfate surface area densities 
and reproduces the observed decrease in ozone following the Mt. Pinatubo eruption. The average loss 
is now greatest in the mid-1 990s instead of around 2000 in simulations without sulfate surface area 
density variability and the extra 5 ppt of Br y . The new simulation also compares better to observations 
during this same period. 

3.2. New ODS Scenario 

The total chlorine prescribed by the mixing ratio boundary conditions is shown in Figure 5a for the A1 2010 
(black curve) and SPARC 2013 (light blue curve) scenarios. The difference in the two scenarios becomes 
apparent after 2040 when the SPARC 2013 total chlorine becomes slightly higher. The SPARC 2013 scenario 
has a maximum difference of 65 parts per trillion by volume (pptv), with 60 pptv more total chlorine than 
the A1 2010 scenario by 2100. Much of this difference is due to the longer lifetime of F1 1. Figure 5b shows the 
annual average concentration of F1 1 in the two simulations in the tropical troposphere. The difference in 
the two scenarios is apparent in F1 1 after about 2015, two decades earlier than it is possible to discern a 
difference in total chlorine. The longer lifetime in the new scenario leads to 14 pptv more F1 1 by 2100; since 
there are three chlorine atoms per F1 1 molecule, this causes 43 pptv more total chlorine. This represents 72% 
of the difference in total chlorine by 2100. 

3.3. Total Column Ozone Response 

Next, we focus on the TCO resulting from three simulations with the latest version (GEOS-5_new) of 
GEOSCCM. The first simulation (P_A1_2010, orange curves) is for the past and includes observed SST and SIC. 



Figure 5. (a) Annual average total chlorine (Cl tot ) in ppbv from the A1 
2010 (black curve) and SPARC_2013 (light blue curve) ODS scenarios. 

(b) Annual average F1 1 concentration (ppbv) in the tropical troposphere 
from the PF_A1_2010 (black curve) and PF_SPARC_201 3 (light blue curve) 
simulations from 1960 to 2099. 
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Figure 6. Annual average quasi-global (60°S-60°N) total column ozone 
(DU) for P_A1_2010 (orange curve), PF_A1_2010 (black curve), and 
PF_SPARC_2013 (light blue curve) simulations. P_A1_2010 simulation is 
from 1960 to 2012, and the PF_A1_2010 and PF_SPARC_201 3 simula- 
tions are from 1960 to 2099. The thick lines are output smoothed to 
remove interannual variability. The black dashed line indicates the 1 980 
value of the smoothed PF_A1_2010 simulation. 


The second simulation (PF_A1_2010, 
black curves) includes both past and 
future but uses modeled SST and SIC 
from an AR5 simulation with CESM1. 
The third simulation (PF_SPARC_2013, 
light blue curve) uses the same SST and 
SIC as PF_A1_2010 but includes the 
new ODS scenario. 

Evolution of quasi-global TCO is similar 
for these three simulations (Figure 6) 
with a few small differences. The 
simulation with observed SST and SIC 
(P_A1_2010) has slightly higher TCO, with 
the largest differences from the other 
simulations between 1990 and 2010. 
Evolutions of TCO are similar for the two 
scenarios with modeled SST and SIC, with 
somewhat less TCO in the second half of 
the 21st century in the PF_SPARC_2013 



Time (yeors) 



Time (years) 

Figure 7. (a) Annual average 10°S-10°N total column ozone (DU) for P_A1_2010 (orange curve), PF_A1_2010 (black curve), and 
PF_SPARC_2013 (light blue curve) simulations. P„A1_2010 simulation is from 1960 to 201 2; the PF_A1_2010 and PF_SPARC_2013 
simulations are from 1 960 to 2099. The thick lines are output smoothed to remove interannual variability. The black 
dashed line indicates the 1980 value of the smoothed PF_A1_2010 simulation. The inset shows TCO and 3°S-3°N 20hPa U-wind. 
(b) Same as Figure 7a but for 25°S-25°N. The inset shows TCO for 25°S-25°N and 20°S-20°N tropical SSTs with reversed scale. 
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Figure 8. Annual average (a) southern polar cap (60-90°S) and (b) northern polar cap (60-90°N) total column ozone (DU) 
for P_A1_2010 (orange curve), PF_A1_2010 (black curve), and PF_SPARC_2013 (light blue curve) simulations. P_A1_2010 
simulation is from 1960 to 2012; the PF_A1_2010 and PF_SPARC_201 3 simulations are from 1960 to 2099. The thick lines are 
output smoothed to remove interannual variability. The black dashed line indicates the 1980 value of the smoothed 
PF_A1_2010 simulation. 


run. These differences are less than 1 DU on average. All three simulations show that the prominent decrease in 
ozone follows the Mount Pinatubo eruption [ Aquila et al., 2013] in 1991. 

Tropical TCO can be significantly impacted by changes in tropical upwelling [Li et al., 2009] as well as by 
changes in ODSs. Lin and Fu [2013] have shown that the Brewer-Dobson Circulation has multiple 
branches including the following: a transition zone between the troposphere and stratosphere, a 
shallow branch between 70 and 30hPa, and a deep branch above 30hPa. Processes such as El Nino- 
Southern Oscillation, or more generally tropical SST anomalies, and the QBO can impact different 
branches of this circulation. This new version of GEOSCCM includes a realistic QBO, which modulates 
tropical upwelling and thereby impacts ozone concentrations. Figure 7a shows the annual average TCO 
over the deep tropics (10°S-10°N) for the three simulations, revealing strong interannual variability. The 
inset shows an expansion of 2070-2090 with the PF_A1_2010 simulations deep tropical TCO (10°S-10°N) 
and the 20hPa U-wind (3°S-3°N), clearly revealing the close relationship of TCO with the QBO in this 
region. The QBO winds induce a secondary meridional circulation, so that deep tropical upwelling 
changes are largely balanced by subtropical changes in the opposite direction. This circulation acts to 
diminish the ozone change due to the QBO when averaged over the broader tropics. Figure 7b shows 
that the interannual variability in TCO due to the QBO is significantly reduced for a wider tropical average 
(25°S-25°N). Over this region tropical SSTs dominate the impact on tropical upwelling variability and 
consequently TCO variability. The inset in Figure 7b shows an expansion of 2070-2090 with the 
PF_A1_2010 simulations broad tropical TCO (25°S-25°N) and tropical SSTs (20°S-20°N) with a reversed 
scale as the two are anticorrelated. This can also be seen in the PF_A1_2010 (black thin curve) and 
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a) Annual Cl y (ppbv) A1 2010 





b) Annual Cl y (ppbv) SPARC 2013 
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Figure 9. The 2050-2080 annual average C\ y (ppbv) for (a) the PF„A1_2010 simulation and (b) the PF_SPARC_201 3 simula- 
tion. The difference in C\ y between these two simulations is shown in (c) pptv and (d) percent. 


PF_SPARC_2013 (light blue thin curve) simulations that show strongly correlated variability forced with 
identical SSTs. During the second half of the 21st century the PF_SPARC_2013 simulation tends to have 
slightly less ozone than the PF_A1_2010 simulation which, like the difference in the quasi-global TCO, is 
less than 1 DU on average. This is again quantitatively consistent with the slightly larger halogen loading 
in the PF_SPARC_2013 simulation. 

Both regions also show slightly larger TCO when observed SSTs are used over modeled ones. This is 
consistent with the slightly reduced tropical upwelling simulated with slightly cooler observed SSTs 
compared to modeled SSTs. In both regions of the tropics TCO does not recover to pre-1980s levels by 
2100 due to the increasing Brewer-Dobson circulation. In GEOSCCM the increase in circulation 
dominates over the impact of cooling stratospheric temperatures in tropical total column ozone 
[Waugh et at., 2009]. For upper stratospheric ozone, cooling dominates [Li et at., 2009], by slowing the 
temperature-dependent ozone loss processes [Haigh and Pyle, 1982]. 

Figure 8a shows the annual SH polar (60-90°S) average TCO evolution for the three simulations. TCO 
evolution from PF-A1_2010 and PF-SPARC_2013 are similar, with no statistically significant differences 
between them. In contrast to the tropical and quasi-global TCO, the P_A1_2010 tends to produce less 
ozone at SFH high latitudes at least prior to 1 995. Figure 8b shows that over the NH polar region (60-90°N) 
the high interannual variability and relatively small difference in ODS scenario do not allow the simulations 
to be distinguishable. An ensemble of simulations would be necessary to see a significant impact from the 
higher halogen loading in the SPARC 2013 scenario at high latitudes. 

3.4. The 2050-2080 Cl y and Ozone 

The differences between the two ODS scenarios are largest during the second half of the 21st century 
when EESC is expected to recover to 1980 levels. On an annual basis ozone is projected to recover 
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a) Annual Ozone (ppmv) A1 2010 



b) Annual Ozone (ppmv) SPARC 2013 
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Figure 10. The 2050-2080 ozone (ppmv) annual average for (a) the PF_A1_2010 simulation and (b) the PF_SPARC„2013 
simulation. The difference in ozone between these two simulations is shown in (c) ppbv and (d) percent. 


somewhat sooner as seen in Figure 8. Figures 9a and 9b show the annual average 2050-2080 Cl y 
concentration (ppbv) for the A1 2010 and SPARC 2013 ODS scenarios, respectively. The difference in the 
two scenarios is shown in Figures 9c and 9d for absolute (pptv) and relative (%) difference, respectively. 
In the upper stratosphere a 60 pptv increase in Cl y is simulated in the PF_SPARC_2013 run, which is 
about a 3-4% difference. This difference is consistent with that seen in l /elders and Daniel [2014]. 

The ozone concentrations (ppmv) for the same time period are shown in Figures 10a and 10b for the A1 2010 
and SPARC 2013 scenarios, respectively. The difference in ozone is less than 50 ppbv (Figure 10c) over much 
of the stratosphere, which is generally less than 2% (Figure lOd). A large ensemble would be required to 
obtain statistically significant differences, especially at high latitudes. 


4. Conclusions 

Improvements to GEOSCCM since Pawson et al. [2008] have on the whole led to significant reductions in 
ozone biases. GEOSCCM now includes an internally generated QBO, which has improved the tropical 
interannual variability. The increase in the maximum solar zenith angle from 90 to 94° as the demarcation of 
day/night chemistry led to a significant reduction in the high-latitude ozone bias. Also, the addition of 5 ppt 
of Br y further reduced the bias especially in the SFH polar region during August and September. 

The GEOSCCM quasi-global (60°S-60°N) total column ozone now compares well to ground-based estimates, 
especially over the last 30 years. Including observed stratospheric sulfate surface area densities, and 
combined with the extra Br y , increases the ozone loss following the Mount Pinatubo eruption and causes the 
largest average ozone loss in the model to occur during the 1990s as observed, instead of near 2000 in the 
previous simulations that lack this sulfate variability. 
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The rate of tropical upwelling increase has decreased moving from the GEOS-4 to GEOS-5 modeling 
framework with identical SSTs and GHG concentrations. The smaller differences between GEOS-5_old_m and 
GEOS-5_new_m are largely due to the different SSTs and GHGs that were used in those runs. 

We used GEOS-5_new to produce simulations with the A1 201 0 and SPARC 2013 ODS scenarios. In the SPARC 
2013 scenario there is 60 pptv increase of total chlorine over the previous scenario by 2100, and about 72% of 
this change comes from F1 1 changes. This difference is most apparent in the second half of the 21st century 
and leads to relatively small differences in the quasi-global ozone, typically less than 1 DU. 

Within the deep tropics (10°S-10°N) the interannual variability is dominated by the QBO. Over the broader 
tropics (25°S-25°N) much of the QBO-induced variability in the deep tropics is balanced by opposing changes 
in the subtropics, leading to relatively smaller net ozone changes. Here tropical SST-induced variability 
dominates. Ozone in the tropics does not return to pre-1980 levels by 2100 due to circulation changes. 
Although the SPARC 2013 scenario increases extratropical Cl y by 3-4% by 2050-2080, its impact on high- 
latitude ozone is much smaller than the effect of interannual dynamical variability. At high latitudes, an 
ensemble of simulations is required to quantify the extra Cl y impact on ozone. 

For 21 st century ozone projections it is important to understand what model biases impact the response to a 
perturbation and which do not. Douglass et al. [2012] demonstrated that background upper stratospheric 
temperatures and NO y values influenced the ozone response to a chlorine perturbation. The changes 
discussed here do not significantly influence these quantities and mainly impact background ozone 
concentrations and its variability. 
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